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NOTATION 


) 


Ij.  Ji.  Kj 


m(zp 


V  T„ 


Si(Zi) 

Xi.  Y..  Z. 


X,  y,  z 


X  ,  y  ,  z 


*0'  yo'  *0 

A.,  B.,  r. 

«.  P.  7 


Amplitude  of  incident  surface  wave 

Radius  of  circle  on  which  spar  axes  are  located 

Radius  of  i^^  spar  as  a  function  of  z^ 

Acceleration  of  gravity 

Moment  of  inertia  of  i^  spar  about  x-,  y-,  z-axes, 
respectively 

2 

Wave  number  of  incident  waves  u  /g 

Mass  of  one  spar  (assumed  equal  for  each  spar) 

Mass  per  unit  length  of  i^  spar 
Nxxmber  of  spars 

Space-fixed  cylindrical  polar  coordinate 
Body-fixed  cylindrical  polar  coordinate 
Defined  in  Equations  [13a],  [13b] 

Cross-sectional  area  of  i^^  spar  -  ra?(zj) 

Components  of  force  on  i^^  spar 
Space-fixed  Cartesian  coordinates 
Body-fixed  Cartesian  coordinates 
Linear  displacements  of  structure 
Components  of  moment  on  i^  spar 

Angular  displacements  of  structure  about  x-,  y-,  z-axes 
Space-fixed  cylindrical  angular  coordinate 
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Body-fixed  cylindrical  angular  coordinate 
Equilibrium  angular  poeition  of  i^^  spar  axis 
Kinematic  viscosity  of  water 
Density  of  water 
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Complete  velocity  potential 
Velocity  potential  of  incident  wave 
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ABSTRACT 


A  theoretical  analysis  is  constructed  for  the  hydro- 
dynamic  forces  acting  on  a  system  of  interconnected  vertical, 
slender,  axisymmetric  bodies  which  are  floating  in  presence 
of  incident  waves.  The  theory  is  based  on  linearized  water 
wave  potential  theory  and  the  use  of  slender  body  techniques. 
The  resulting  expressions  for  the  hydrodynamic  forces  are 
used  to  predict  the  motions  of  such  a  system.  The  effects 
of  viscous  damping  are  also  estimated. 


INTRODUCTION 

A  spar  raft  as  defined  here  consists  of  several  long  thin  bodies  of 
revolution  rigidly  interconnected  so  that  they  will  float  vertically  in  the 
water  and  support  a  platform  or  submerged  weight.  When  regular  waves 
are  incident  on  such  a  structure,  it  will  generally  oscillate  in  six  degrees 
of  freedom.  The  purpose  of  this  report  is  to  provide  an  approximative 
method  for  calculating  such  motions. 

The  assumptions  are;  (1)  that  the  spars  are  identical,  (2)  that  their 
interconnections  are  made  in  such  a  way  that  the  mass  and  the  hydrody¬ 
namic  effects  of  the  connecting  members  may  be  neglected,  and  (3)  that 
the  individual  spars  are  far  enough  apart  for  their  hydrodynamic  interac¬ 
tions  to  be  neglected.  The  motions  of  a  single  spar  buoy  have  been  treated 
by  Newman;^  here  his  method  is  extended  to  the  case  of  N  spars  arranged 
in  a  circle.  In  addition  to  including  the  hydrodynamic  and  inertial  forces 
on  several  spars,  it  is  necessary  to  extend  Newman's  analysis  to  allow 
for  all  six  degrees  of  freedom.  (In  his  problem,  only  three  degrees  of 
freedom  involved  nontrivial  results.) 

The  basic  assumptions  of  Newman's  analysis  are  used  here.  In 
particular,  it  is  assumed  that  the  wave  amplitudes  and  body  motions  are 
small  enough  that  linearized  free  surface  theory  may  be  applied  and  that 


^References  are  listed  on  page  30. 


the  spar  radii  are  small  enough  compared  to  wavelength  and  spar  sepa¬ 
rations  that  slender  body  theory  may  be  used.  Equations  of  motion  are 
derived  on  these  bases.  These  equations  predict  motions  which  are  un¬ 
damped;  thus  they  are  valid  only  for  frequencies  which  are  not  near  the 
resonance  frequencies. 

Near  the  resonance  frequencies,  it  is  necessary  to  consider  the 
damping  due  to  wave  generation,  and  this  report  shows  that  forces  of  high¬ 
er  order  in  terms  of  spar  radii  must  be  included.  The  leading  damping 
forces  are  found,  thus  providing  equations  of  motion  valid  near  resonance. 

In  addition,  this  report  indicates  that  viscous  forces  depend  linearly 
on  the  velocities  for  axial  motions,  and  these  forces  are  found  explicitly. 


GEOMETRY  AND  COORDINATES 

It  is  convenient  to  define  several  coordinate  systems.  With  the 
structure  floating  at  rest,  we  place  the  origin  of  a  space-fixed  reference 
frame  at  the  undisturbed  free  surface  over  the  center  of  gravity  of  the 
structure.  Let  the  Cartesian  coordinates  of  a  point  in  this  system  be 
(x,y,  z),  with  the  z-axis  directed  upwards.  In  this  same  system  we  de¬ 
fine  cylindrical  coordinates  (R,  6,  z); 

r2  =  +  y2  ;  e  =  tan"^  - 

'  X 


or 

X  =  R  cos  0  ;  y  =  R  sin  0 

z  is  here  the  same  as  the  Cartesian  coordinate  z. 

Let  the  undisturbed  axis  of  the  i^^  spar  be  located  at  R  =  aQ,  0  =  0^. 
We  define  another  set  of  space-fixed  coordinates,  (x^,  y^,  zj),  with  origin 
at  R  =  aQ,  0  =  0p  z  =  0.  Let  the  cylindrical  coordinates  of  a  point  in  this 
system  be  (R^,  z^),  with  the  latter  having  the  same  orientation  as  the 

previous  cylindrical  system. 

In  the  undisturbed  condition,  the  surface  of  the  i^^  spar  will  be 
specified  by  the  equation; 
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t'Vi 

Figure  1  -  The  Spar  in  Its  Equilibrium  Position 

R.  =  a.(zj)  ;  zj  >  -  Hi 

These  quantities  are  all  shown  in  Figure  1. 

Finally,  we  introduce  primed  coordinate  systems  which  correspond 
to  each  of  the  systems  just  mentioned,  but  which  are  fixed  in  the  body. 
When  the  body  is  in  its  equilibrium  position,  the  primed  and  unprimed 
systems  coincide. 

The  linear  displacement  of  the  raft  will  be  described  by  three  dis¬ 
placement  variables,  Xq,  yQ,  Zq  ;  the  angular  displacement  by  a,  P,  V» 
which  are  the  positive  rotations,  respectively,  about  the  x-,  y~>  z-axes. 
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Since  the  motions  are  assumed  to  be  small  enough  that  squares  and  prod¬ 
ucts  of  these  variables  are  negligible^  the  location  of  the  raft  is  complete¬ 
ly  specified,  and  the  three  angular  displacement  variables  can  be  treated 
as  the  components  of  a  vector. 

Let  be  the  position  vector  of  a  point  fixed  in  the  raft,  where 
(x^,y^,z^).  In  terms  of  the  space-fixed  coordinates,  r  =  (x, y,z), 
we  have,  to  first  order  in  small  qua'htities: 

£  =  £o  +  2^£' 


where 


lo  =  <*o*  yo'  *o) 

a  =  (or,  p,  Y  ) 

In  terms  of  components,  this  equation  is  equivalent  to 

x=  x'+XQ+Pz'-Yy'’ 

y  =  y'+  yo+  vx'-az"  [i'] 

z  =  z'  +  Zq  +  Qfy'-px' 

To  first  order  also,  it  follows  that 

x'  =  x-XQ-Pz  +  yy 

y'  =  y-yo-y*+  [^"1 

z'  =  z-ZQ-ory+Px 
The  unprimed  coordinates  are  related  by 

X  =  Sq  cos  0^  +  X£ 

y  =  ap  sin  0.  +  y^  [2] 
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and  the  same  equations  hold  if  x,  y,  z  and  y*i  z-  are  all  primed. 

It  is  assumed  that  there  are  incident  waves  which  are  described  by 
the  velocity  potential 

(x,y,z,t)  =  e*'*  cos(kx- wt)  [3] 

where  A  is  the  amplitude  of  surface  wave, 

k  is  2r/ wavelength  =  wave  number,  and 
w  is  the  circular  frequency. 

In  the  definition  of  the  coordinates  given  above,  the  orientation  of  the  x- 
and  y-axes  was  not  specified,  except  for  the  orientation  of  the  plane  which 
they  defined.  Now  we  specify  that  the  x-axis  points  in  the  direction  of 
propagation  of  the  surface  waves  and  the  y-axis  completes  the  right-hand 
system. 


THE  VELOCITY  POTENTIAL 

The  surface  of  the  i^^  spar  can  be  specified  by  the  equation: 
0  =  F,(x,-,  y,'  z/) 

^2  ,  ^2  t  ww 

=  X.  +  y.  -  a.  (Z.  ) 

=  [xi  -  Xq  -  Pzj  +  Y(ao  sin  +  y^)]^ 

2 

+  [yi  "70“  Y(ao  cos  6^  +  xj)  +  az^] 

-  aj  [z.  -  Zq  “  ®i  yi^  +  ®i 

fh 

The  boundary  condition  on  the  i''  spar  is  then 
8F 

-~  +  (Vi«-V.)F.  =  0  on  F.  =  0 


5 


where  Vj  indicates  the  gradient  in  the  (x^.y^.z.)  system,  and  #s*(x.,y.,  z.,t) 
is  the  velocity  potential  (viz.,  plus  a  potential  due  to  the  presence  and 
motion  of  the  structure).  After  some  simplification,  we  find  that  the 
boundary  condition  is 


89  _  ^  89 

9Ri  dzi 


[xq  +  P(z.  +  a.  ap  -  Sf  ao  sin  ej  cos  X. 

+  [yQ  -  i(z.  +  a^ap  +  Y  Sq  cos  ej  sin  X. 


[4] 


-  [zq  +  aQ(of  sin  0^  -  P  cos  6^)]  on  =  a^ 

where  a.  =  -  .  Second  and  higher  order  terms  in  the  motion  variables 

1  dz^ 

have  been  consistently  dropped. 

Let  ^  substitute  this  relation  into  the  last  equation 

to  obtain  a  boundary  condition  on  .  From  Equations  [2]  and  [3]  we 
note  that 


-  (jjA  kz  .. 

*0  =  e  cos  (kx  -  wt) 

^  k 

[5] 

=  cos  [k(aQ  cos  cos  \^)  -  wt] 

k 


It  follows  that 

w  Ae^^i  ^  +  ka^  a^  +  f  a?)  cos 

+  I  k^  cos  sin  (kag  cos  0.  -  wt) 

-  [(a/*  +  j  ka^)  +  7  ka.  cos  ZX^]  cos  (kaQ  cos  -  wt) 

+  [61 

where  the  omitted  terms  are  of  third  and  higher  order  in  terms  of  the 
spar  radius  a*  and  its  derivative  or  of  second  and  higher  order  in 


\0Ri  i  ®*i/Fj=0 
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the  motion  variables.  Neglecting  now  second-order  terms  in  a.  and  a/*, 
we  find  the  condition  on 


/8«j  f 


kz*  ^ 

+  (a)  A  e  ^  sin  (kaQ  cos  B*  -  wt)  r  cos 

+  ^Yq  -  iz.  +  y  ap  cos  }  sin 

+  ka  ju)  cos  (kaQ  cos  0^  —  cjt)J  cos  2X 

-  ^Zq  a-'  +  aQ  a^'  (dr  sin  -  P  cos  8^) 


[7] 


1  kz  *  ^ 

-  (a/'  +  j  ka. )  w  Ae  ^  cos  (kaQ  cos  0^  -  wt)  > 


We  note  that  the  boundary  condition  is  now  applied  at  the  surface  of  the 
undisturbed  spar,  and  the  right-hand  side  is  evaluated  in  the  space-fixed 
(xp  y^,  z^)  coordinate  system. 

If  this  i^^  spar  were  located  alone  in  an  infinite  fluid  with  the  above 
boundary  condition  valid  for  -  H  <  z^  <  0,  the  solutions  for  by  slender 
body  theory  would  be 


[  zgaj'  +  aga^' (a  sin  0^  -  P  cos  0j) 

-  a)Ae^^  (a^''  +  j  ka^)  cos  (kaQ  cos  0^  -  wt)] 

+  |a?[xQ  +  p^-Y3LQ  sin  0^  +  w  A  e^^  sin  (kaQ  cos  0^  -  cot)] 


+  i  a?  [yQ  -  +  Y  aQ  cos  0.]  — 

-  i  a^  [wkAe^^  cos  (kag  cos  0^  -  wt)]  +  (Zj  -  1)^]  ^  d; 

BxfJ 
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,  dai(U 

where  ai  -  a'(^)  and  a-  ~  "  ■  ■■  in  the  integrand.  We  adapt  this  type  of 

*  ^  *  d4  1 

solution  to  the  present  case  as  follows:  (1)  [Rj  +  (z*  -  4/  1  ^  poten¬ 

tial  for  a  source  located  at  (0,  0,  C)  in  an  infinite  fluid. 

We  replace  this  source  potential  by  another  source  potential  which, 
in  addition,  satisfies  the  free  surface  condition.  (Z)  We  now  impose  the 
condition  that  the  radius  a^  is  much  smaller  than  the  distance  between 
spars;  i.e.,  a^/aQ  «  1 .  Then  the  potential  obtained  by  satisfying  the 
boundary  condition  on  the  i^^  spar  will  produce  negligible  fluid  velocities 
at  the  other  spars,  and  the  total  potential  can  be  expressed  as  a  sum  of 
potentials,  each  satisfying  the  conditions  on  one  spar.  The  resulting  total 
potential  is 


^(x,  y,  z,  t)  =  e^^  cos  (kx  -  wt) 

k 

^  Q  f 

+  ^  J  aj  [zQa/  +  aQa/  (or  sin  0.  -  P  cos  0^) 

n=l 

-  o)  Ae^^  (a/  +  j  ka^)  cos  (kaQ  cos  0^  -  wt)] 


+  1  a?  [xq  +  p  ^  “  Y  ^0  0^  +  w  A  sin  (kaQ  cos  0^  -  wt)] 


+  i  [VO  "  +  V  ^0  cos  0j]  ^ 


-  I  a^  [wkAe^^  cos  (kaQ  cos  0j  -  wt)] 

|[Rf  +  (z.  -  ^  Jq(  vRi)dv|  d^ 

N  0 

+  JTuk  ^a.  [zga.'  +  aQa." (a  sin  0^  -  P  cos  0^) 


[8] 
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[8] 

.  .  continued 

+  +  j  ka^)  sin  (kaQ  cos  0^  -  wt)] 


+  a 


?  [xq  +  ~  Y  ®i  ®i  ” 


+  af  [Yq-  “C  +  Y^o  cos 


+  ^  a^[k  A  sin  (k  aQ  cos  0^  -  ut)] 


In  accordance  with  the  assumptions  of  slender  body  theory,  we  eval¬ 
uate  these  terms  as  0  and  identify  the  values  so  obtained  with  the 

potential  on  the  body  surface  =  a^.  Because  a^/aQ  «  1,  the  value  of  ^ 
to  lowest  order  on  the  i^^  spar  depends  only  on  the  first  term  above  and 
one  term  in  the  first  sum.  By  the  same  approximation  procedure,  we  find 
that  all  terms  in  the  second  sum  contribute  amounts  of  higher  order  in 
terms  of  a^.  Later  we  shall  reconsider  the  second  sum  when  we  calculate 
damping  forces. 


FIRST-ORDER  FORCES  AND  MOMENTS 
The  pressure  is  obtained  from  Bernoulli's  equation  in  linearized 

form 


P  =  -  P  ^  -  Pgz  [9] 

Thus  it  will  be  necessary  to  evaluate  on  each  spar  and  to  integrate, 
in  an  appropriate  way,  the  result  over  all  spars  to  find  the  forces  and 
moments. 

Using  slender  body  approximations  again,  we  find  that  on  the  i^^ 

spar 
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#  =  6*^*^  COS  (kaQ  cos  0^  -  wt)  -  2«  a^  cos  X-  sin  (kaQ  cos6^  -  ut) 

-  [zpa.'  +  aQa.'lo  sin  0.  -  P  cos  0^ 

-  uAe^^Ma."  +  j  ka.)  cos  (kaQ  cos  0.  -  wt)]  a^  log  a. 

% 

-  [xq  +  Pz.  -  vaQ  sin  0j  a.  cos  X.  -  [y^  -  az-  +  ya^j  cos  0.  ]  a.  sin  X. 

+0(a?)  [10] 

Also,  we  note  that  on  the  surface  of  the  spar 

z  =  z/"  +  Zq  +  Qr(aQ  sin  0^  +  sin  Xp  -  P(aQ  cos  0.  +  a^  cos  X." ) 

Thus  the  pressure  on  the  i^^  spar  is 

f  kz'^ 

p  =  -pJgAe  1  sin  (kaQ  cos  0,  -  wt) 

kz  ^ 

+  2gkAe  ^a.  cos  cos  (k  aQ  cos  0^  -  wt) 

..  [in 

-  [  Zq  a/  +  aQ  a/ ( a  sin  0.  -  P  cos  0^  ) 

kz  ^ 

-  gkAe  1  (a/  +  j  ka.)  sin  (kaQ  cos  0^  -  wt)]  a.  log  a. 

-  [xq  +  pz.'-  yaQ  sin  0.]a.  cos  X.'  -  [y^  -  az(  +  y  a.^  cos  0.]a.  sin  X." 

+  g  [z/"  +  Zq  +  a(aQ  sin  0.  +  a^  sin  X/)  -  P(aQ  cos  0^  +  a^  cos  Xj')]^ 

Let  the  force  be  resolved  along  the  space-fixed  axes  which  corre¬ 
spond  to  the  coordinates  (x,y,z).  In  particular,  designate  the  components 
of  hydrodynamic  force  on  the  i^^  spar  by  Xj,  Y^,  Z^.  Likewise,  let  the  com 
ponents  of  hydrodynamic  moment  on  the  i^^  spar  be  denoted  by  A^,  B. , 
which  correspond  to  the  rotations  a,  p,  V*  Note  specifically  that  the 
moments  are  taken  with  respect  to  the  space-fixed  axes  at  the  center  of 
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the  whole  structure.  If  we  let  n  be  a  unit  normal  vector  out  of  the  fluid, 
then 


Xi=  J 

Si 

p  co8(n, x)dS 

[12a] 

V,=  J 

p  co8(n,y)dS 

[12b] 

J 

p  cos(n, z)dS 

[12c] 

J 

p  [y  cos  (n,  z)  ~  z  cos  (n,  y)]  dS 

[I2d] 

Bi=  J 

p  [z  cos  (n,  x)  -  X  cos  (n,  z)]  dS 

[I2e] 

=  1 

f_ 

p  [x  cos  (n,  y)  -  y  cos  (n,  x)]  dS 

[I2f] 

Si 


The  integrals  are  taken  over  the  instantaneous  surface  of  the  i^  spar. 
Here  cos(n,  x)  is  the  cosine  of  the  angle  between  n  and  the  x-axis,  etc. 
We  find  readily  that,  to  first  order  in  small  quantities, 

cos(n,  x)  =  -  cos  X/  +  y  sin  X/  +  Pa/' 

cos(n,  y)  =  -  -y  cos  X/'  -  sin  X/'  -  ora/' 

cos(n,  z)  =  P  cos  X-"  -  a  sin  X-"  +  a/' 

For  abbreviation,  we  also  define  two  sets  of  integrals 

s„  =  /"  S(.p(.,-)“dZi' 

•'-H 

T„  =  r°  eW  S(zp(2.')"dz.-' 


[13a] 

[13b] 
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2  tK 

where  8(21")  =  ra^  (z*')  s  cross-sectional  area  of  spar  at  z^. 

The  combination  of  these  formulas  and  definitions  with  the  previous 
pressure  results  yields  for  the  force  and  moment  components 

Xj  =  2  pgkA  cos(kaQ  cosB.  -  wt)TQ  -  p(xq  -  y  Sq  sine.)SQ  -  pPSj  [  14a] 

Yj  =  -  p(yo  +  yao  cos  0i)So  +  paSi  [14b] 

Z-  =  pgkA  sin(kaQ  cos  0.  -  wt)TQ  +  pgSQ 

[14c] 

-p  g  [  A  sin(kag  cos  -  wt)  +  Zq  +  aaQ  8in0j  -  P  aQ  cos  0j]  8(0) 

A.  =  pgkaQ  sin  0|  A  8in(kaQ  cos  0.  -  cjt)TQ  [  14d] 

+  Pg(yO  ^0  +  Y  ao  cos  +  p(yQ  +  y  cos  0.  -  g»)Sj  -  paS2 

-  p  g  aQ  sin  [  Zq  +  sin  0^  -  (3  aQ  cos  0^  +  A  sin(k  aQ  cos  0*  -  wt)]  S(  0) 

-  p  g  kaQ  cos  0^  A  sin(kaQ  cos  0^  -  wt)  Tq  +  2  pg  k  A  cos  (kaQ  cos  0^-  wt ) 

-  Pg(^0  ■*“  ^0  •  VaQ  sin  0.)Sq  -  p(xq  -  y  aQ  sinO.  +  gP)S2  -  p  pS2 

+  p  g  aQ  cos  0^  [zq  +  aaQ  sin  0-  -  p  aQ  cos  0^  +  A  sin  (kaQ  cos  0^  -  wt)]  S(0) 

[14e] 

r.  =  -  2  pgkaQ  sin  0-  A  cos(kaQ  cos  0-  -  wt)TQ 
+  p  [ag  sin  0.  Xq  -  ag  cos  0J  yo  -  a^v]  Sq 

[14f] 

+  p[aQ  sin  0^  P  +  aQ  cos  B^a]  S| 

S(0)  is  the  cross-sectional  area  at  z^  =  0  when  the  whole  system  is  at 
rest  and  in  equilibrium* 


t 


SIMPLIFIED  INTERPRETATION  OF  FIRST-ORDER 
FORCES  AND  MOMENTS 


These  expressions  for  the  forces  and  moments  can  be  viewed  from 
a  simple  point  of  view.  Consider,  for  example,  the  x-component  of  force, 
Equation  [I4a],  which  when  written  out  becomes 


Xj  =  p  r  S(z.')  /z  gkAe^*i  cos  (kan  cos  6-  -  wt) 

•'-H 

-  Xq  +  Y  aQ  sin  -  z/^^dz/  [14a'] 


From  Equation  [5]  we  see  that 

®0  kz  ' 

-  =  gkAe  ^  cos  (kan  cos  0.  -  wt) 

0t8x  ®  0  1  ' 

fVi 

on  the  equilibrium  position  of  the  i  spar.  Thus  the  first  term  in  the 
bracket  in  Equation  [I4a‘^]  is  just  twice  the  local  acceleration  that  the 
water  would  have  at  the  mean  position  of  the  spar  axis  if  the  spar  were 
not  present.  The  terms,  -  xq  +  Y  ^  •  give  the  negative  ac¬ 

celeration  (in  the  x-direction)  of  the  point  on  the  spar  axis.  The  quantity 
pS(Zj^)  is  the  added  mass  per  unit  length  of  a  cylinder  accelerating  nor¬ 
mally  to  its  axis.  Thus  the  x-component  of  force  is  the  integral  over  the 
mean  spar  length  of 

(added  mass  per  unit  length)  times  (2  times  local 
water  acceleration  on  spar  axis  due  to  incident  wave 
alone  minus  acceleration  of  point  on  axis  of  spar). 

It  may  appear  strange  that  the  water  particle  acceleration  is  doubled 
in  this  formula.  However,  the  cause  is  seen  on  examination  of  Equations 
[8]  and  [10] .  In  the  latter  equation,  the  terms  containing  the  factor 
(cos  X.-)  give  rise  to  x-components  of  force.  Here  the  term  due  to  the 
incident  waves  (the  second  term)  is  already  doubled.  Half  of  this  con¬ 
tribution  comes  directly  from  the  first  term  of  Equation  [8]  (i.e.,  direct¬ 
ly  from  the  incident  wave  potential)  and  half  comes  from  the  term 
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a?  wA  e 


sin{kao  cot  0.  -  cot)  ^  [R?  +  (z.  -  ^4 


The  latter  is  effectively  a  diffraction  potential;  it  is  part  of  the  singularity 
potential  which  offsets  the  normal  velocity  component  of  the  incident  wave 
on  the  spar.  These  results  can  also  be  regarded  as  a  special  case  of  a 
general  body,  accelerating  in  a  time-varying  (but  spatially  constant)  in¬ 
finite  field  of  fluid.  It  follows  from  consideration  of  the  forces,  both  in 
the  fixed  and  moving  coordinate  systems,^  that  the  hydrodynamic  force 
on  the  body  is  the  added  mass  times  the  relative  acceleration  plus  the 
displaced  mass  of  fluid  times  the  spatial  acceleration  of  the  (undisturbed) 
fluid.  For  a  circular  cylindrical  section,  the  added  mass  and  the  dis¬ 
placed  mass  are  equal,  and  the  above  relation  for  the  x- component  of  the 
force  on  the  spar  follows  immediately.  • 

The  z-component  of  force,  Equation  [14c],  consists  of  three  parts: 

(a)  [  pg  Sq]  ; 

(b)  [“  pg(zQ  +  aaQ  sin  0^  -  paQ  cos  0^)  S(0)]  ; 

(c)  pg  A  sin  (kaQ  cos  0^  -  tijt)  [k  Tq  -  S(0)]  . 


Part  (a)  is  just  the  hydrostatic  force.  Part  (b)  is  the  decrease  in  buoyancy 
which  occurs  when  the  spar  is  raised  an  amount  (zq  +  aaQ  sin  0^  -  P  aQCos  0^). 
Part  (c)  is  the  integral  over  the  undisturbed  spar  surface  of  the  vertical 
pressure  force  due  to  the  incident  wave  alone.  This  is  easily  seen  by 
noting  that 


kTQ-S(O) 


and,  thus,  that  Part  (c)  is  equal  to 


pgAe 


kzjj' 


sin  (ksQ  cos  0j[  ~  (a)  t)] 


Pq  cos  (n,  z)  dS 
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to  first  order,  where  Pq  is  the  first  term  of  Equation  [11]  • 

The  moments  are  obtained  by  calculating  the  force  per  unit  length 
along  each  spar,  multiplying  by  the  appropriate  lever  arm,  and  integrating 
along  the  lengths  of  the  spars.  It  should  be  noted  again  that  the  moments 
are  calculated  with  respect  to  the  space-fixed  axes.  Thus  a  point  located 
at  on  the  i^^-axis  has  space-fixed  coordinates  (see  Equation  [I"*])* 

X  =  Sq  cos  0.  +  Xq  +  p  z/  -  \  Sq  sin  0^  ; 

y  =  aQ  sin  +  yo  +  V  cos  0^  -  az-' ; 

z  =  z*^  +  Zq  +  oaQ  sin  0*  -  Psq  cos  0*  . 


FIRST-ORDER  EQUATIONS  OF  MOTION 

Let  M  be  the  mass  of  a  spar.  (The  N  spars  are  assumed  to  be  iden¬ 
tical.  )  Let  I^,  be  the  moments  of  inertia  of  the  i^^  spar  about  the 

X-,  y-,  z-axes,  respectively.  Moreover,  let  Mq,  Iq,  Jq,  and  Kq  denote 
the  mass  and  moments  of  inertia  for  any  additional  superstructure,  and 
assume  (0,  0,Zq)  to  be  the  center  of  gravity  thereof.  Then  the  equations 
of  motion  are 

N 

(Mq  +  NM)iQ  =  Xi  [15a] 

i=l 

N 

(Mo  +  NM)yo  = 

i=l 
N 

(Mq  +  NM)2q  =  27  Z.  -  g(NM  +  Mq) 
i=l 

N  N  N 

-  Mq  g  zq  +  g  ^  ^  -  g  ^  f  m{z()y  dz(  [15d] 

i=0  i=l  i=l  ^ 


[15b] 

[15c] 
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N  N 

2* 

i=0 

N 

£' 

i=0  i=l 


N 


i=l  i=l 

miz^x  dz." 

[15e] 

N 

[15f] 

m(z*^)  is  the  mass  per  unit  length  of  the  spar  itself.  The  integral  is  taken 
over  the  length  L  of  the  spar.  This  length  generally  extends  from  =  -H 
to  some  value  of  z^  greater  than  zero,  x  and  y  are  the  distances  from  the 
fixed  reference  frame  to  a  point  on  the  axis  of  the  i^^  spar. 

The  moment  of  inertia  about  the  x-axis  is 

li  =  [  m(zj'){y2  +  z2)dz{’ 

•'l 

where  y  =  aQ  sin  +  yo  +  V  S-o  -  crz^'^  +  ...»  and 

z  =  Zj*  +  Zq  +  or  aQ  sin  0^  -  P  aQ  cos  +  ...» 

the  omitted  terms  being  of  higher  order  in  the  small  motion  variables. 
Since  is  multiplied  by  5,  we  need  keep  only  the  zero-order  terms  in 
1^.  Clearly  then. 


to  the  required  order  in  small  quantities;  that  is,  has  the  same  value  as 
in  the  equilibrium  position.  Similarly, 

=  I  m(z(){x^  +  z^)dz/ 

=  [  miZi")  [ao  cos^Gi  +  z(^]  dz{ 
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I 


X 


Ki  ■  mCz/Xx^  +  y^)d*/' 

=  f  m(z/')aQdzj  =  aQ  M 
•'L 

The  integral  terms  in  Equations  [I5d]  and  [I5e]  can  also  be  written 
explicitly  by  expanding  x  and  y  in  the  integrands.  Thus 

J  m(z.')y  dz.'  =  J  m (z.')  [ag  sin  0.  +  y^  +  y  ag  cos  0.  -  or  z/]  dz." 

Li  Li 

=  (aQ  sin  +  Yq  +  Y  ^0  cos  0j)M  -  «  J  m(z.')z.‘'dz.‘' 

Xj 

J  m(z.'')  X  dz.'  =  m(z.')  [aQ  cos  0j  +  Xq  +  p  z."  -  y  aQ  sin  0.]  dz/ 

=  (aQ  cos  0j  +  Xq  -  y  aQ  sin  0j)M  +  P  J  m(zpz.'' dz." 


We  note  that  the  final  integral  terms  here  would  have  vanished  if  z."  had 
been  measured  from  the  spar  center  of  gravity. 

Now  let  us  write  the  equations  in  full.  Equation  [15a]  becomes 


N 

(Mq  +  NM)xq  =  2  pgkATQ  cos  (kaQ  cos  0.  -  ut)  -  pNxqSq 

i=l 


N 


[15a'] 


-  pNp  Si  +  paQy  Sq  2’  sin  6^ 


Clearly,  Mq  +  NM  =  pSqN,  since  at  rest  the  buoyancy  of  the  total  raft 
equals  its  weight.  Also,  we  now  impose  the  condition  that  the  spars  have 
a  regular  angular  spacing.  Thus 


Zir 


ej  =  0J+  ;  i  =  1.2, ...,N 
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If  N  s  1,  then  Eq  -  0,  and  the  laet  term  in  Equation  [l5a‘']  vanishes. 
If  N  >  1. 


N-1  2iri 

^  sin  0i  =  ^  sin  (ej  +  — j  = 
i=l  i=0 


•  U  ■  N-l  2ir\  ,  /N  2r\ 


.  0 


•‘“S 


and  again  the  last  term  in  [iSa"*]  vanishes. 

Therefore,  the  equation  of  motion  for  xq  is 

N 


2(Mq  +  NM)xq  +  pNSj  M2  pgkATo  2"  cos  (kaQ  cos  0.  -  wt)  [l6a] 

i=l 


By  similar  arguments,  the  equation  of  motion  for  yq,  Equation 
[l5b],  becomes 


2(Mq  +  NM)yQ  -  pNSj  a  =  0 


[16b] 


The  equation  for  Zfx,  written  out,  is 


N 


(Mq  +  NM)(zq  +  g)  =  pgkATo  2  sin  (k  Uq  cos  0.  -  wt)  +  N  p  g  Sq 


i=l 

N 


-  pg  AS(0)  sin  (kaQ  cos  0^  -  wt)  -  pgNzoS(O) 


i=l 


N 


-  pgaaoS(O)  2  sin  0j  +  pgpaoS(O)  2  cos  0J 


i=l 


N 

r 

i=l 


Again  we  note  that  Mq  +  NM  =  pSgN,  which  enables  us  to  eliminate  the 
gravity  term  on  the  left.  Also, 


2 

4-1 


COS  0^ 


=  0 


N  >1 
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and  aQ  =  0  for  N  =  1 .  The  equation  becomes 

N 

(MQ*fNM)zQ  =-N  p  g  Zq  S(0)  -  p  g  A  [-kTQ  +  S(0)]  sin  (kap  cos  0^  -  <ot) 

i=l 

[16c] 

The  first  term  on  the  right-hand  side  is  just  the  change  in  buoyancy  which 
accompanies  a  vertical  displacement  of  the  raft.  The  remaining  terms 
correspond  to  the  vertical  force  obtained  by  integrating  the  dynamic  pres¬ 
sure  due  to  the  incident  wave  over  the  surface  of  the  spar.  This  is  the 
Froude-Krylov  hypothesis;  To  a  first  approximation,  the  presence  of  the 
body  does  not  distort  the  incident  wave  or  the  pressure  associated  with  it. 

The  or- equation  is 


.  N  N  . 

«  1^4  ^  +  Iq  +  ^  dz.'l 


N 

=  -  pgaoA[S(0)-  kTo]  sin  0^  sin  (kaQ  cos  0^  -  a>t) 

i=:l 


+  NpSjy-Q-NpgofSj  -  NpaS2 


N 

-  p  g  aQ  S(0)  [a  sin^  0^  -  P  sin  0.  cos  0.  ] 

i  =  l 


+ 


m(zpz/  dz/  +  MqZq 


I 


I 


If  N  >  2, 


i: 

{si 


sin^  0j 


|N 
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Thus,  for  N  >  2, 


2*  sin  cos  =  0 
i=l 

o  "^jNM  ag  +  N  p  S2  +  Iq  +  N  J  m  (zplz-')^  dz-'^  +  a  p  g  Sj  [I6d] 

+  iNpgaQS(O)  -  Mggzg  -  Ng  J  m(z.')z.'dz/^ 

N 

-NpSiyo  =  -  pgaoA[S(0)-kTo]  sin  0^  sin  (kaQ  cos  0^  -  wt) 


Similarly,  for  N  >  2,  the  p-equation  becomes 


p  ^iNMaJ  +  NpS2  +  Jg  +  N  dz.'|  +  p^NpgSj 

+  iNpgaQS(O)  -  Ng  m  { zp  z.' dz.' 

N 

+  N  pSj  Xg  =  pgagA  [S(0)  -  kTg]  ^^cos  0^  sin  (kaQ  cos  0-  -  wt) 


IN 

+  2  pgkATj  cos  (kaQ  cos  0^  -  ojt) 


[16e] 


Under  the  same  assumptions,  we  obtain  for  the  last  equation 

N 

^Kq  +  2NMaQ^  =  -  2pgkaQATQ  sin0^  cos  (kaQ  cos  0.  -  wt)  [16f] 


In  the  case  of  N  =  1,  the  above  equations  reduce  to  Newman's  equa¬ 
tion  for  a  single  spar.  If  N  =  2,  these  equations  do  not  hold.  However, 
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the  special  conditions  that  follow  from  N  =  2  can  easily  be  applied  here  to 
obtain  simple  equations.  The  case  is  not  considered  sufficiently  interest¬ 
ing  to  warrant  writing  out  the  equations  here* 

The  heaving  motion  can  be  obtained  immediately^  if  desired,  since 
the  ZQ-equation  contains  no  coupling  terms.  In  addition,  the  equation  for 
rotation  about  the  z-axis  is  not  coupled  to  the  other  equations.  However, 
the  Yq  and  a  motions  are  coupled;  also  the  Xq  and  P  motions.  Similarly, 
the  couplings  are  simple  enough  for  these  equations  to  be  solved  directly. 

We  should  note  that  there  are  two  resonance  frequencies.  In  heave, 
there  is  resonance  when 


2  pgS(O) 


M  +  Mq/N 


[17a] 


In  either  of  the  coupled  motions  there  is  resonance  when 

ZMgj^pSj  +  ApaQS(O)  -  MqZq/N  m(z/')z/  dz/j 


0) 


ZM^lMaJ  +  pS2  +{ko/n}+  Jl 


[17b] 


Since  the  equations  contain  no  damping  forces,  infinite  response  ampli¬ 
tudes  are  predicted  when  resonance  occurs.  Of  course,  this  is  meaning¬ 
less  in  the  linearized  model  and  so  the  above  equations  ([I7a]  and  [17b]) 
can  be  valid  only  in  frequency  ranges,  not  including  neighborhoods  of  the 
two  exceptional  frequencies.  When  such  neighborhoods  are  excluded  from 
consideration,  the  predictions  should  be  fairly  accurate  if  the  small  ampli¬ 
tude  and  slenderness  restrictions  are  observed,  since  damping  forces  are 
of  higher  order  than  the  forces  considered.  Near  the  resonance  frequen¬ 
cies,  however,  the  damping  forces  are  important,  even  if  small.  This 
problem  is  considered  in  the  next  section. 
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EQUATIONS  OF  DAMPED  MOTION 


In  the  previous  section,  we  considered  only  buoyancy  and  accelera¬ 
tion  forces  on  the  spars  because,  generally,  these  were  the  forces  of 
lowest  order  in  the  small  parameter  a^(z).  At  resonance,  these  forces 
cancel  each  other,  thus  they  are  no  longer  the  lowest  order  forces.  We 
must  re-examine  the  previous  analysis  and  include  terms  of  a  higher 
order  in  a^(z)  to  obtain  equations  of  motion  which  have  meaning  at  and 
near  resonance. 

The  boundary  condition,  Equation  [7],  was  valid  only  to  first  order 
in  a^(z).  If  we  now  include  second-order  quantities  (in  a^tz))  in  the  bound¬ 
ary  condition.  Equation  [7],  and  add  the  necessary  corresponding  terms  in 
the  potential  function,  Equation  [8],  we  simply  obtain  more  terms  in  the 
acceleration  and  buoyancy  forces.  Since  these  terms  are  much  smaller 
than  those  already  considered,  they  can  alter  the  response  only  slightly, 
principally  by  changing  the  resonance  frequencies  somewhat.  They  still 
contribute  no  damping  forces. 

Nevertheless,  the  desired  damping  forces  can  be  obtained  from  the 
second  summation  in  the  potential  function,  Equation  [8].  These  terms 
were  discarded  earlier  because  they  contributed  forces  of  higher  order 
than  those  being  considered.  It  is  easily  seen,  however,  that  these  terms 
do  lead  to  damping  forces,  which  will  be  the  lowest  order  forces  at  reso¬ 
nance.  We  also  see  that  the  terms  in  this  second  sum  which  involve  the  in¬ 
cident  wave  amplitude  A  do  not  contribute  damping  effects.  They  simply 
affect  the  driving  force,  again  by  an  amount  of  Ixigher  order. 

So  now  we  consider  the  potential 


N  Q 

=  jTwk  j  -^a.  [zQaj‘'+  aQa.'(a  sin  Oj  -  p  cos  0j)] 

i=l 


+  a?  [xq  +  p  ^  -  Y  ao  sin  0.  ] 


+  a?[yQ-a^  + Vaocos  0.]  Jo(kR.)]  d; 


. 
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(See  the  second  sum  in  Equation  [S], )  Here  and  are  functions  of  4* 

th  8®* 

To  calculate  the  associated  force  on  the  i  spar,  we  must  evaluate  - 

at 

for  R.  =  aj(Zi). 

For  small  values  of  (kR>),  the  Bessel  function  in  can  be  approx¬ 
imated  by  the  beginning  of  its  Taylor  series,  that  is, 

Jo(kRj)  =  1  -  i(kR.)2  +  ... 

Similarly, 

^  Jo(kRi)  =  -ik^R.  cos  X.  +  ... 

—  Jo(kRi)  =  -  ik^R.  sinX.  +  .  .  . 

dyi 

Keeping  a  one -term  approximation  in  each  case,  we  find 


=  iwk[zQ  +  &Q{a  sin  0.  -  p  cos  6.)]  [S(0)  -kT^] 

Ri=ai 

-  i<ok^a.(z.)e^^i  |[(xq-  sin  e.)TQ  +  pTjcos  X. 

+  [(Yq  +  Y  ag  cos  Sj)  Tq  -  orTj  ]  sin  XjJ. 


The  pressure  due  to  this  potential  is,  when  evaluated  on  R^  =  a^. 


Ri=ai 


-  p 


Ri=ai 


=  -  2pwk[zQ  +  aQ(«  sin  0.  -  P  cos  0. )]  [S(0)  -  kTg] 

+  i  p«k^a.(z.)e'^*^  "  V  ag  sin  ei)Tg  +  pTj]  cos  X. 

+  [(Yq  ^*0 
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In  the  expressions  for  and  several  terms  which  are  of  higher 
order  in  small  variables  have  been  omitted. 

The  forces  and  moments  due  to  this  pressure  distribution  are  calcu¬ 
lated  from  Equations  [12a]  through  fl2f],  with  asterisks  inserted  where 
appropriate.  The  results  are  as  follows: 

xf  =  -ipa,k^{(io-Vao  8ine.)T2  +  pToTj 

Yf  =  -  ip,,k^^(yO  +  Yao  cos  0. )  tJ  -  i  Tq 

Z*  =  -  I  pwk  ^Zq+  aQ(<i  sin  6.  -  P  cos  e.)^[S(0)  -  kT^]^ 

A*  =  -  A  pwkaQ  sin  0^  -[zq  +  aQ(Q'  sin  0.  -  P  cos  0.)^  [S(0)  -  kTjj]^ 

+  i  pwk^{{yO  +  VaQ  cos  6i)TQTj  - 

B*  =  I  pwkaQ  cos  0j  >[zq  +  aQ(<i  sin  0^  -  P  cos  0^)^  [S(0)  -  kTp]^ 

-  ipu)k^{(iQ-  Y^o  sin  0i)TQTj  +  pT^j- 

There  is  no  need  for  a  damping  moment  F*,  since  the  Y^i^otion  has  no 
resonance  in  any  case. 

N 

The  modified  equations  of  motion  are  obtained  by  adding  ^  xf, 

i=l 

etc.,  to  the  right-hand  sides  of  the  previous  equations,  [I5a]  through 
[l5f],  or  alternatively,  to  [l6a]  through  [l6f].  After  simplification, 
the  equations  become >  for  N  >  2, 


N|[2{M  +  Mo/N)xQ  +  ipa)k^TjiQ]  +  [pS^P  +  i  p cok^  Tq p]^ 


N 

=  ZpgkATQ  cos  (kaQ  cos  0j  -  wt) 
i=l 


[18a] 
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N  |[2(M  +  MQ/N)yo  +  ipwk^xJyQ]  -  [pS^a  +  |  pwk^  Tq  Tj  «]J>  »  0  [18b] 


N^(M  +  MQ/N)zQ  +  ipcok[S(0)-kTQ]2  Zq  +  pgS(0)zQ]- 

N 

=  -  pgA[S(0)-  kTQ]  z  sin  (kaQ  cos  6.  -  wt) 


N  {(iMaQ  +  pS2  +  mz^dz]  a  pukaj  [S(0)  -  kT^j]^ 

+  j  pwk^Tjja  +  [pgSj  +  j  pgaQS(O)  -  gMQZQ/N-gJ  mzdz] 

Lj 

-  [pSilVo-  [iP<*>k3ToTj]yo} 
N 

=  -  pgaQA[S(0)-  kXo]  2"  sin  sin  (kaQ  cos  6.  -  ut) 


N  {[iMaJ  +  pS2  +  Kq/N+I  mz^dzjp  +/i  p«kaJ[S(0)-  kXg]^ 

Lj  ' 

+  |pwk^Xjjp  +  [  pg  Sj  +  i  pgaQS(O)  -  gMQZQ/N-gJ*  mzdz]p 

Li 

+  [pSi]io+  [ipwk^ToXjlxQ^ 

N 

=  pgaoA[S(0)- kXol  2  cos  6^  sin  (kaQ  cos  6^  -  ut) 

i=l 

N 

+  2  pgkAXj  2  COS  (kaQ  cos  6.  -  a>t) 


[18c] 


[I8d] 


[I8e] 
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N 

(Kq  +  NMaQ)  Y  =  -  pgkaQATQ  sin  0^  cos  (kaQ  cos  0.  -  iot)  [I8f  ] 

i=l 

Of  course  these  equations  can  be  solved  very  easily  by  substitution 
of  Xq  =  Cj  sin(a)t  +  6|),  yQ  =  C2  8in(wt  +  62),  etc.  The  results  will  not 
be  written  because  no  additional  perspicuity  seems  to  follow. 


VISCOUS  DAMPING 

All  damping  forces  introduced  so  far  correspond  to  the  energy  lost 
through  radiation  of  surface  waves.  In  addition,  energy  will  be  lost  through 
the  mechanism  of  viscosity.  The  viscous  damping  forces,  in  general,  will 
be  of  second  order  in  the  motion  variables.  As  an  example,  suppose  that  a 
right  circular  cylinder  translates  in  a  direction  perpendicular  to  its  axis. 
The  viscous  drag  is  proportional  to  the  velocity  squared,  and  so  is  negli¬ 
gible  by  the  standards  already  assumed. 

If  the  cylinder  has  an  axial  motion,  however,  the  viscous  force  will 
be  linear  in  the  velocity.  For  example,  if  a  right  circular  cylinder  of 
radius  a  has  an  axial  velocity  Re  {We^^^}  =Wq  cos(cjt  +  €  ),  then  we  can 
see  easily  from  elementary  fluid  mechanics  that  the  velocity  anywhere  in 
the  fluid  is 


where  K0(x)  is  the  modified  Bessel  function  of  argument  x  and  order 
zero,  and  the  second  expression  gives  Kq(x)  in  the  Kelvin  representation 
of  its  real  and  imaginary  parts.  The  only  velocity  component  is  that  which 
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is  parallel  to  the  cylinder  axis,  and  this  component  depends  in  space  only 
on  the  distance  r  from  the  axis. 

A  tangential  (axial)  stress  on  the  cylinder  surface  is  given  by 


9w 

9r 


r=a 


Wob 


cos 


^wt  +  €  +  p  -  Y  + 


where  b,  c,  p,  y  are  real  numbers  defined  by 


be^^  =  keri  n/^  a  +  ikeii  s/^a 
1  V  Ay 

c  e^^  =  ker  a  +  i  kei  a 


Thus  there  is  an  axial  force  component  per  unit  length  on  the  cylinder: 


F  =  - 


2  TT  p  a  nTcjv  WQb 


cos  ^wt  +  €  +  p  -  Y  + 


Tj 


2  TT  pa  n/wv  WQb 


^cos(a)t  +  c)  cos  ^P  -  Y  + 
(wt  +  «)  sin  ^p  -  Y,+ 


—  Sin 


The  first  term  in  brackets  in  the  last  equation  gives  the  drag  force. 

For  the  slender  bodies  considered  in  this  report,  it  is  consistent  to 
calculate  the  viscous  force  in  a  stripwise  manner.  That  is,  at  each  cross 
section  of  the  i^^  spar,  where  the  radius  of  the  section  is  a^(zp,  we  con¬ 
sider  that  particular  section  to  be  part  of  an  infinitely  long  right  circular 
cylinder  translating  axially,  calculate  the  axial  viscous  force  per  unit 
length,  and  integrate  such  results  over  the  length  of  the  spar.  To  first 
order  in  the  motion  variables,  the  axial  velocity  of  the  i^  spar  is 

Zq  +  ^0(0  sin  ~  ®i) 


27 


Thus  we  add  to  Equations  [14c],  [I4d],  and  [I4e],  respectively,  the  addi¬ 
tional  force  and  moments 


2  r  p  ^ycov 


a^(z.)b.(z.) 

Ci(Zi) 


COS 


*i)  -  Y(Zi)  + 


(zq  +  aQOf  sin  0^  -  aQ^  cos  0j)dZ| 


+ 


ai(Zi)b.(z.) 

UCi(Zj) 


sin 


P(Zi)  -  Y(zj)  + 


Bf*  = 


ao  sin  e.  Zf* 

-  aQ  COS  0.  Z]*'* 


+  ano 


sin  0^  -  aQ^  COS  0^)dz.^ 


If  these  additional  quantities  are  inserted  into  the  equations  of  mo¬ 
tion,  we  can  still  obtain  solutions  by  the  same  method  used  previously. 

Although  the  viscous  forces  thus  obtained  are  linear  in  the  veloci¬ 
ties,  they  do  not  fit  properly  into  the  perturbation  scheme  in  terms  of 
small  radius.  The  modified  Bessel  functions  encountered  have  singulari¬ 
ties  when  the  argument  approaches  zero.  In  fact. 


b  ^  1 

a  a  log  a 


as  a  -►  0. 


Thus,  as  the  slenderness  of  the  spars  is  accentuated,  the  viscous  forces 
increase.  This  is  in  contrast  to  the  potential  flow  forces,  which  become 
smaller  and  smaller.  No  general  conclusion  can  be  drawn  concerning  the 
relative  importance  of  the  viscous  and  nonviscous  damping  forces  as  far 
as  dependence  on  radius  is  concerned.  Calculations  should  be  made  in 
individual  cases. 
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DISCUSSION 


Equations  of  motion  have  now  been  obtained»  based  on  the  assumption 
that  wave  amplitude »  body  motions,  and  spar  radii  are  small  quantities.  If 
damping  is  neglected,  there  are  two  resonance  frequencies  (as  given  by 
Equations  [I7a]  and  [17b])  at  which  infinite  amplitudes  of  motion  will  oc¬ 
cur.  Of  course,  there  does  exist  damping  which  prevents  the  responses 
from  actually  being  infinite.  That  part  of  the  damping  due  to  generation  of 
outgoing  waves  is  small  (of  second  order  in  terms  of  spar  radii)  and  so 
is  of  importance  only  near  resonance.  Generally  there  is  also  a  viscous 
damping  which  may  be  of  importance  throughout  the  interesting  range  of 
frequencies,  or  possibly  only  near  resonance,  or  perhaps  not  at  all.  This 
damping  is  associated  only  with  axial  velocities  of  the  spars,  but  its  ef¬ 
fects  appear  in  each  of  the  modes  of  motion  for  which  there  is  a  resonance. 
The  relative  importance  of  viscous  damping  can  be  determined  in  individ¬ 
ual  cases  only  by  actually  carrying  out  solutions  of  the  equations  of  motion. 

If  it  is  desired  to  minimize  the  motions  of  the  structure  over  a  range 
of  practical  wave  lengths,  the  most  effective  procedure  is  to  attempt  to  re¬ 
move  the  natural  frequencies  from  the  desired  range.  If  this  is  not  possi¬ 
ble,  the  natural  frequencies  should  be  chosen  as  frequencies  for  which  the 
incident  wave  amplitudes  are  smallest.  In  practice,  this  will  generally  be 
equivalent  to  reducing  the  natural  frequencies  as  low  as  possible. 

From  Equation  [  I7a]  we  see  that  the  heave  natural  frequency  is 
proportional  to  the  radius  of  the  spars  at  the  undisturbed  waterline.  (S(0) 
is  the  cross-sectional  area  at  the  equilibrium  waterline.)  Assuming  that 
the  total  mass,  Mq  +  NM,  is  approximately  fixed,  we  have  no  other  param¬ 
eter  to  adjust;  thus  we  would  minimize  S(0)  as  far  as  possible. 

In  the  case  of  the  other  resonance  frequency,  we  see  from  Equation 
[l7b]  that,  apparently,  there  are  several  parameters  available  for  adjust¬ 
ment.  We  note  that  if  the  center  of  buoyancy  and  center  of  gravity  coin¬ 
cide,  then 


pSi  -  MqZq- 


I  in(z)  z  dz  =  0 

•'l 
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By  raising  the  center  of  gravity  and/or  lowering  the  center  of  buoyancy* 
we  can  make  this  difference  negative  and  thereby  possibly  decrease  this 
natural  frequency.  However*  in  general*  the  exact  amount  by  which  such 
adjustments  will  affect  the  natural  frequency  is  not  certain*  since  the  de¬ 
nominator  in  Equation  [17b]  will  also  be  affected.  In  any  case,  reductions 
in  S(0)  (for  lowering  the  heave  frequency)  will  also  lower  the  pitch-surge 
frequency. 
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